Cellular buckling in I-beams 



Cellular buckling from mode interaction in 
I-beams under uniform bending 

By M. Ahmer Wadee*and Leroy Gardner 

^vj I Department of Civil & Environmental Engineering, 

Imperial College London, London SW7 2AZ, UK 

^. 

Abstract 

»j> ' Beams made from thin-walled elements, whilst very efficient in terms of the structural 

strength and stiffness to weight ratios, can be susceptible to highly complex instability phe- 
r^ . nomena. A nonlinear analytical formulation based on variational principles for the ubiquitous 

I-beam with thin flanges under uniform bending is presented. The resulting system of differ- 
ential and integral equations are solved using numerical continuation techniques such that the 
response far into the post-buckling range can be portrayed. The interaction between global 
lateral-torsional buckling of the beam and local buckling of the flange plate is found to oblige 
the buckling deformation to localize initially at the beam midspan with subsequent cellular 
»vj ' buckling (snaking) being predicted theoretically for the first time. Solutions from the model 

Qs ' compare very favourably with a series of classic experiments and some newly conducted tests 

r~^ . which also exhibit the predicted sequence of localized followed by cellular buckling. 



> 



^ 



^ ■ 1 Introduction 

Beams are possibly the most common type of structural component, but when made from thin 
metallic plate elements they are well known to suffer from a variety of elastic instability phenomena. 
There has been a vast amount of research into the buckling of thin- walled structural components 
^\ ' with much insight gained ilj |2] i3j . In the current work, the well-known problem of a beam made 

from a linear elastic material with an open and doubly-symmetric cross-section - an "I-beam" - 
under uniform bending about the strong axis is studied in detail. Under this type of loading, long 
beams are primarily susceptible to a global mode of instability namely lateral-torsional buckling 
(LTB), where, as the name suggests, the beam deflects sideways and twists once a threshold critical 
moment is reached 4 . However, when the individual plate elements of the beam cross-section, 
namely the flanges and the web, are relatively thin or slender, elastic local buckling of these may 
also occur; if this happens in combination with global instability, then the resulting behaviour 
is usually far more unstable than when the modes are triggered individually. Other structural 
components, usually those under axial compression rather than uniform bending, such as I-section 
struts with thin plate flanges [S] , sandwich struts [5] , stringer-stiffened and corrugated plates [71 |H] 
and built-up, compound or reticulated columns fS] are all known to suffer from so-called interactive 
buckling phenomena where the global and local modes of instability combine nonlinearly. 

Experimental evidence suggests that the destabilization from the mode interaction of LTB and 
flange local buckling is severe |10[ 111] , the response being highly sensitive to geometric imperfec- 
tions particularly when the critical loads coincide |T21[I3]. Nevertheless, apart from the afore- 
mentioned work where some successful numerical modelling (both finite strip and finite element) 
and qualitative modelling using rigid links and spring elements were presented, the formulation 
of a mathematical model accounting for the interactive behaviour has not been forthcoming. The 
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current work presents the development of a variational model that accounts for the mode interac- 
tion between global LTB and local buckling of a flange such that the perfect elastic post-buckling 
response of the beam can be evaluated. A system of nonlinear ordinary differential equations sub- 
ject to integral constraints is derived, which are solved using the numerical continuation package 
Auto [H]. It is indeed found that the system is highly unstable when interactive buckling is 
triggered; snap-backs in the response showing sequential destabilization and restabilization and 
a progressive spreading of the initial localized buckling mode are also revealed. This latter type 
of response has become known in the literature as cellular buckling |15| or snaking |16] and it is 
shown to appear naturally in the numerical results of the current model. As far as the authors 
are aware, this is the first time this phenomenon has been found in beams undergoing LTB and 
local buckling simultaneously. Similar behaviour has been discovered in various other mechanical 
systems such as in the post-buckling of cylindrical shells [17l [18] and the sequential folding of 
geological layers [HI [20] . 

Experimental results generated for the current study and from the literature are compared 
with the results from the presented model; highly encouraging quantitative results emerge both in 
terms of the mechanical destabilization exhibited and the nature of the post-buckling deformation 
with tangible evidence of cellular buckling being found. This demonstrates that the fundamental 
physics of this system is captured by the analytical approach. A brief discussion is presented on 
whether other structural components made from thin-walled elements may also exhibit cellular 
buckling when local and global modes of instability interact. Conclusions are then drawn. 



2 Variational formulation 

2.1 Classical critical moment derivation via energy 

Consider a uniform simply-supported (pinned) doubly-symmetric I-beam of length L made from 
an isotropic and homogeneous linear elastic material with Young's modulus E and Poisson's ratio 
v. The overall beam height and flange widths are h and b respectively with the web and flange 
thicknesses being t^, and t/ respectively. The beam is under bending about the strong axis with a 
uniform moment M, as shown in Figure [Ija). The second moments of area about the strong and 
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(a) Beam geometry: elevation (left) and cross-section (right) 




(b) LTB displacements 



Figure 1: (a) An I-section beam under uniform strong axis moment M. (b) Lateral displacements 
and torsional rotations due to LTB. A right-hand Cartesian coordinate system is employed. 



weak axes are defined as I^ and ly respectively. When LTB occurs, as noted in the literature [4], 
the strong axis bending moment and corresponding displacements only couple with the lateral dis- 
placements and torsional rotations at higher orders. For this reason and for the sake of simplicity, 
the displacements arising from strong axis bending are presently neglected, which is a common 
assumption; in future work, however, these coupling effects may be incorporated. Therefore, only 
two separate lateral displacements are defined for determining the respective positions of the local 
centroids of the web and the flanges: Us and u^,. The lateral displacements of the top (ut) and 
bottom flanges (ub) respectively are thus (see Figure [TJb)): 



Ut 



Ub 



(1) 
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Moreover, a torsional angle of magnitude (p to the vertical is introduced and the applied moment M 
is transformed thus into components of strong and weak axis moments, M^. and My respectively: 



M, 



Mcos{-(j>) « Af, My = Msin(-(/)) « -M0, 



(2) 



both expressions being written to the leading order. Figure[2ja) shows a plan view of the top flange 
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(a) Axial stresses due to bending 



<^hx 





O-0 
(b) Bending moments under LTB 



Figure 2: After LTB: (a) top flange stress components from strong (abx) and weak axis (cFby) 
bending moments; (b) induced bending moment components during LTB, M^ and My, which 
force one of the flange outstands to become more heavily compressed and hence vulnerable to 
local buckling. 

of the beam and the stress distribution components from the strong axis and the induced weak axis 
bending moments. Therefore, the maximum strong axis compressive stress occurs when y ~ h/2 
and coupling this with the co-existing compressive component of the weak axis stress, the most 
vulnerable flange outstand is defined, as shown in Figure [2l^b). Since the important component of 
bending under LTB is about the weak axis, the values of the relevant second moment of area for the 
flange (//) can be approximated to Iy/2 of the whole section, where ly = b^tf/6, by assuming that 
the overall contribution of the web is very small, which is true for Lbeams of practical dimensions. 
The strain energy stored in the beam due to bending (Ub) is therefore: 



Ub = ^Elf 



dz = -EL, 



uf dz + I EL. 



(3) 



where primes denote differentiation with respect to the axial coordinate z, Uw — h(j)/2 and /„ 
lyK^ jA. The strain energy stored from uniform torsion {Ut) is: 



Ut = 7: 



Td0== -GJ 
9 / 9 



J,'2 



dz. 



(4) 



where T is the torque, G is the material shear modulus and J is the St Venant torsion constant 
defined for an Lsection as: 

'"-'+3 I n^ o+ .^+3 1 (•g\ 



J 



[2btj + ih-2tf)ti 



The work done by the external moment MO is given by the induced weak axis moment My 
multiplied by the average end rotation from bending about the weak axis; this is given by the 
following expression [21] : 



Me 



and so the total potential energy is thus: 



V^Ub + Ur- MQ = 



—EI,,u 



J/"s 



M(j)u'^ dz 



-EIu 
2 



-GJc^'^ - Mcjju'^ 



dz. 



(6) 



(7) 



To find the critical moment Ma-, the calculus of variations could be used to derive the governing 
differential equation. However, since the solution of the buckling mode is known to be sinusoidal 
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[4] , the same result can be achieved by using a two degree of freedom Rayleigh-Ritz formulation 
with the following trial functions: 



qsLsin 



L' 



-90! 



L' 



(8) 



where Qs and q^ are generalized coordinates defining the amplitudes of Ug and respectively. Note 
that owing to the coordinate system used, the negative sign in front of q^ is necessary to ensure 
that (?0 > when qg > 0. Substituting the expressions for Us and from Equation ([5]) gives: 
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GJtt^ 
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-90' 
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dz 



(9) 



where A = h/2L. The formulation is in small deflections (linear) and so only a critical equilibrium 
analysis is possible at this stage. The advantage of using the Rayleigh-Ritz method becomes more 
apparent in the interactive buckling model below. Performing the integration and then assembling 
the Hessian matrix V^j at the critical point C, gives the following condition: 



vH- 



•^ ss 

vF 



Vi 
v^ 

y A.A. 



-AU,{-K/Lf GJ{-K/Lf + EI^{'k/LY 



where the individual elements of the matrix V~ are thus 



= 0, 



(10) 



^ ah 



d'^v 



(11) 



Solving equation ([T0| gives the classical expression for the critical moment Mcr that triggers LTB 
for a beam with a doubly-symmetric cross-section under uniform bending ^ : 



M„^-yi^;G7yi + ^f^ 



(12) 



2.2 Interactive buckling model 

From the previous section, it has been shown that as the displacements and rotations from LTB 
grow the applied moment M can be expressed as a component about the strong axis (M^;) and an 
induced component about the weak axis (My) at any point along z. As a result, the vulnerable 
outstand of the flange, as identified in Figure [21 may therefore buckle locally as a plate. The 
critical stress of plate buckling for a uniaxially and uniformly compressed rectangular plate, with 
one long edge pinned and the other free, is given by the well known formula [4j: 



atcai = 0.43 



-k'D 



h)i, 



(l-i) 



where, in the current case, 6/ is the width of the vulnerable flange outstand and is given by 
(6 — iu,)/2 and D is the fiexural rigidity of the fiange plate that is equal to £'i'^./[12(l — i^^)]. This 
addresses the case for the flange buckling locally before any LTB occurs. 

It was shown in [6^ that the intrinsic assumptions in Euler-BernouUi bending theory were in- 
sufficient to model any interaction between global and local buckling modes. The allowance of 
shear strains to develop within the individual elements, however small, being key to the formu- 
lation. Figure [3] shows a useful way that shear can be introduced; the displacement of the top 
and bottom flanges being decomposed into separate "sway" and "tilt" modes [22 after the global 
instability (LTB) has been triggered. Each original generalized coordinate q^ and q^ is defined 
as a "sway" and has a corresponding "tilt" component, with associated generalized coordinates 
qt and q^- respectively. This is akin to a Timoshenko beam formulation |23| : when qg ^ qt and 
14> ¥^ It, shear strains are developed and allow the potential for modelling simultaneous LTB and 
local buckling. 
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(a) Top flange 



"/t = -{* - grAJTricos — 
0„, ss — -— = Iq, - g^A)7rcos-- 



(b) Bottom flange 



Figure 3: Sway (upper row) and tilt (lower row) modes of LTB in both flanges. 



Previous work on this type of interactive buckling has included experimental work combined 
with effective width theory [lOj, some phenomenological modelling using rigid links and springs 
along with experiments |llj , some numerical work using a finite strip formulation |,f 2; .f 3 , and a 
finite element formulation [SI]. To model this analytically, however, two displacement functions 
to account for the extra in-plane displacement u and out-of-plane displacement w (Figure |4]) need 
to be defined. Since the strain from the weak axis moment is linear in x and that the boundary 
condition for the line of the web is pinned out of the plane of the flange, the following linear 
distribution in x can be assumed for u and w: 



2x 2x 

u{x,z) = — —u(z), w{x,z) — — —w(z). 



(14) 



It is worth noting that restricting both interactive buckling displacement functions u and w to 
the vulnerable part of the compression flange confines the current model to cases where LTB 
occurs before local buckling. In the cases where local buckling occurs first, the system would be 
expected to trigger the global mode rapidly afterwards [25] and the current model can be used 
to indicate the deformation levels where the interaction would occur (as seen later). However, to 
obtain an accurate linear eigenvalue solution for pure local buckling in the current framework, at 
least another set of in-plane and out-of-plane displacement functions, replicating the role of u and 
w respectively, would need to be defined for the non-vulnerable part of the compression flange; 
this addition is left for future work. 



2.2.1 Local bending energy 

Experimental evidence, presented in [24^, suggests that during interactive buckling the significant 
local out-of-plane displacements within the flanges are confined to the vulnerable outstand. The 
component of additional strain energy stored in bending Um is hence given by: 
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(15) 



Article submitted to The Royal Society 



T^ Paper 



M. A. Wadee and L. Gardnt 



\ 





Web line 


\ 


/ 


IX 


/ 


\\ 


6/2 / 


' "\ 


/ 


1 \ 


/ 



Free edge 




(a) u(x, z) 




2x 
w{x, z) = — 7-w{z) 



(b) w(x, z) 



Figure 4: Displacement functions u and w defined for the vulnerable part of the top flange: 
X = [-6/2, 0] and y = [h/2 - tf, h/2]. 



Substituting w into Uti , the following expression is obtained: 



D 



L ru2 



-w"^ + 4(1 - v)w'^ 
6 



dz. 



(16) 



2.2.2 Flange energy from axial and shear strains 

Since the flanges are assumed to behave in the manner of Timoshenko beams, the bending of 
the flanges when LTB occurs introduce both axial and shear strains, e and 7 respectively. The 
vulnerable part of the top flange, where x — [—6/2,0] and y =- [h/2 — tf,h/2], also has the 
possibility of local buckling; von Karman plate theory gives a standard expression for the axial 
strain e in the z-direction that accounts for both LTB and local buckling terms, thus: 



eti 



duff du 1 f dw 

A h 

dz 



dz 



dz 



[qt + qrX)— sin — 



TXZ 2x 



2x2 



■w ■ 



(17) 



For the part of the top flange that is not vulnerable to local buckling, where x ~ [0, 6/2] and 
y = [h/2 — tf, h/2], and the bottom flange, where y = [—h/2, —h/2 + tf], the following respective 
axial strain expressions are obtained: 



£t2 

£b 



-5— = (<?* + 9rA)— — sm— , 
dz L L 



dutb 
dz 



{qt - qrX)—sm—. 



(18) 
(19) 



The standard strain energy expression is then integrated over the volume of the flanges: 
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dz, 



(20) 



where ip = h/ L is a beam aspect ratio parameter. It is worth noting that the transverse dis- 
placement and the strain in the x-direction are omitted from the current formulation. This is a 
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simplification derived from [7], wfiere tliese components were found to have a negligible effect on 
the post-buckling stiffness of a uniaxially compressed long plate with one longitudinal edge being 
pinned and the other being free. 

In terms of the shear strain 7 in the xz plane within the top flange, where y — [h/2 — tf, h/2], 
von Karman plate theory gives a standard expression, which needs to account for both LTB and 
local terms for the vulnerable outstand x = [--b/2, 0]: 

7ti = Sst — 0u + local buckling terms 

dust dutt du dw dw ,, , , . , ttz 2 Ax . (21) 

and purely LTB terms for the non- vulnerable part of the top flange, x = [0, 6/2], and the bottom 
flange, where y = [—h/2, —h/2 + tf], respectively: 

dUst dutt TTZ 

7*2 = Ost - Ott = -^ h -J— = [{qs - qt) + {q<i> - grjAJ ttcos — , (22) 

dUsb dutb ,. . , . , TTZ 

lb = Osb - &tb = -^ 1- ^^ = [{qs - qt) - (90 - 9r)AJ TTCOS — . (23) 

oz ox Lj 

After obtaining the shear strains, the standard strain energy expression needs to be integrated 
over the volume of the flanges, thus: 
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h/2 j-O fh/2 pb/2 f-h/2+tf fb/2 
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dz 



^ l2 ( "^ + o'^^^'^ + ^^^' ) - His - qt) + (90 - '7r)A] ^ COS Y (2" + ww') 



+ [{qs - qtf + (90 - qr)^X'^] 271^ COS^ ^ I dz. 



(24) 



2.2.3 Work done contribution 



An additional contribution from the vulnerable flange's local in-plane displacement function u 
needs to be included in the work done term. Figure [5] shows that the compression flange has a 
distribution of in-plane displacement which is assumed to have an average linear distribution in x. 
Including this as an average end-rotation, angles ao and a^ are obtained as shown in the diagram. 
Since the end rotation angles can be expressed in terms of the local in-plane displacement function 
M, the expression for the local contribution to the work done is: 

MQi = A4(ao - aL) = -— / u' dz. (25) 



2/1.70 

2.2.4 Total potential energy 

The expression for the total potential energy of the interactive buckling model system can be 
written as a sum of the individual terms of the strain energies minus the work done terms from 
H2.ll and ij2.21 with [/,„ and Ug replacing Ub to account for the change in the bending theory 
assumptions: 

v = u,n + UT + Ubi + Us- M{e + ei). (26) 

This new energy function V replaces equation ([T]) and is written in terms of non-dimensionalized 
variables that replace the original ones, thus: 

~ 2 ^2 ^2 

u = —u, w — yuj, z — —z, (27) 

Ij Ij Ij 
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Average flange end-displacement 
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Outstand not vulnerable to local buckling 



m(OK 
Beam elevation: (x = 0) 



Outstand vulnerable to local buckling 
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Bottom flange 




Figure 5: Contribution from local buckling in the work done term, where h is the overall cross- 
section height. 

giving the full expression for V ^ where primes henceforth represent derivatives with respect to z: 
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(28) 



2.2.5 Linear eigenvalue analysis 

With u and w being zero, along with their derivatives, before any local buckling occurs, the Hessian 
matrix V^, now including terms associated with the "tilt" generalized coordinates, can still be 
used to find Mcr the critical moment for LTB. The Hessian matrix is thus: 



Vf 



with the individual terms being: 
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(29) 



V = Vrc^ = ^GLbtjn'X', Vsr = Vrs = Vtr = Vrt = Vt^ = V<pt = 0, 



(30) 
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which when substituted into V^ with the singular condition at the critical point C, where M 
M„, gives: 



Mr. 



TT EhGJ 



'1 



1 7r2 EL, 



1 + s V (l + s)L2 GJ 



(31) 



This new expression for Men replaces equation (J12p and is used subsequently. The term, s = 
ETr'^tp'^/12G, accounts for the non-zero shear distortion of both flanges, which tends to zero if G 
or L become large; this is an entirely logical result reflecting the difference between Timoshenko 
and Euler-Bernoulli beam theories [23]. However with s > 0, the above expression gives values 
that are marginally below those given by the classical critical moment given in equation (I12p . 



2.2.6 Equilibrium equations 

The total potential energy V with the rescaled variables can be written thus: 

y— £{111" , w , w , u , u; z) dz . 
Jo 



(32) 



Equilibrium equations are found where V is stationary and this is established from setting the 
first variation of V (or SV) to zero, where: 
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(33) 
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Su > dz. 



The term in the integral has to vanish for all Sw and (5w, which gives two coupled nonlinear 
ordinary differential equations: 



6(1 - ly) .„ 3GLHf , 



3ELHj 
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= 0, 



(34) 



3 

u" + -w'w" + 



— \ Tn/' cos — \{qs -qt + A(90 - 9r)) - s{qt + gr A)] - iu+ -ww' ] [ = 0, (35) 

subject to the following boundary conditions which arise from minimizing the terms outside the 
integral in equation p3p : 



;(0) = w"{Q) = wi2) = w"i2) = 0, 



'(0) + ^*'^(0) + 



3A/ 



= 0, 



Ebhtf 
u(l) = w'(l)=w"'(l) = 0, 



(36) 

(37) 
(38) 



where equation p6p refers to pinned boundaries, equation p7p refers to the end strain condition, 
and equation p8|) refers to reflective symmetry of w and antisymmetry of u about the midspan 
respectively. The symmetry conditions are particularly pertinent when LTB occurs simultaneously 
with or before local buckling owing to the sinusoidal distribution of (p forcing the maximum bending 
stress to be located at midspan. 
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Other equilibrium equations can be obtained by minimizing the energy with respect to the 
generahzed coordinates Qs, g^, Qt and qr- 
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3 Physical experiments 

3.1 Specimens and procedure 

A series of physical experiments were conducted on I-beams fabricated by spot- welding thin- walled 
channel sections, made from cold-formed steel, back-to-back. The key material properties were 
measured to be thus: Young's modulus E = 205 kN/mm^, Poisson's ratio ly = 0.3 and the yield 
stress CTy = 290 N/mm^. The channel sections were 75 x 43 x 2 mm in terms of depth, width and 
thickness respectively - see Figure IHl^a). The actual geometry (including corner radii etc.) was 



2 mm 43 m^ 



2 mm, 



80 mm 




38 mm 
(a) Actual (measured) section geometry 




71 mm 



38 mm 
(b) Idealized section geometry 



Figure 6: Cross-section of steel beams tested in the experimental programme, (a) Actual section 
constructed from two channel sections; (b) idealized section geometry derived from experiments. 



converted into an idealized I-section comprising only flat plate elements based on the mean slope 
of the linear regions of the measured load versus maximum bending displacement curves from the 
tests. The dimensions of h and b were hence adjusted slightly such that a meaningful comparison 
with the theory could be made - see Figure [6l^b). Each beam, the idealized properties of which 
are given in Table [TJ had an overall length of 4 metres, was tested under four-point bending with 
a specified buckling length Le given in Table [2] that was controlled by an adjustable pair of lateral 
restraints (see Figure [7] in §3.2p . The critical mode was determined by comparing the strong axis 
bending stress when M 



Mc,., thus: 



Mcr2/n 



'LTB 



2/x ' 



(43) 
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Geometric Property 


Value 


Cross 


-Section Property 


Value 


b 


80 mm 




4 


4.81 X lO'^ mm"* 


h 


71 mm 




ly 


1.71 X 10^ mm^ 


tf 


2 mm 




I^ 


2.16 X 10* mm^ 


^w 


4 mm 




J 


1.86 X 10^ mm^ 



Table 1: Idealized section properties of the experimental samples. 



Test 


Effective length Lg (mm) 


^Local/'''LTB 


ffC (N/mm2) 


Critical mode 


1 


3200 


1.19 


186 


LTB 


2 


3200 


1.19 


186 


LTB 


3 


3000 


1.10 


202 


LTB 


4 


2750 


0.98 


221 


Local (marginally) 


5 


2500 


0.87 


221 


Local 


6 


2250 


0.75 


221 


Local 



Table 2: Buckling lengths for each beam test which reflect the symmetric position of the lateral 

s; these lengths were chosen such that cr^^^^i/ 
triggering both modes approximately simultaneously. 



restraints; these lengths were chosen such that o'Locai/'^'LTB ~ [0.75, 1.20] with one beam (Test 4) 



with the critical stress of plate buckling being given by equation (fT3 



3.2 Testing rig 

A schematic of the experimental setup along with an idealized representation is shown in Figure [71 
The total applied load 2P was split into two point loads each of P applied at a distance xl from 
the end supports. Hence, from simple statics, the uniform moment M between the two loading 
points was Pxl. The loading was displacement controlled; it was applied with a hand-operated 
hydraulic jack in conjunction with a gravity load simulator, a mechanism that adjusted the position 
of the load application relative to the deflecting beam such that the applied load remained vertical. 
Since the jack was hand operated, the displacement was applied in short controlled increments 
but it did mean that dynamic loads were inevitable to a small extent. At midspan, the vertical 
displacement of the beam and lateral displacements of the flanges were measured using linear 
variable differential transformers (LVDTs); the locations of which are presented in Figure |51 

The large displacements and dynamic behaviour of unstable post-buckling, even with displace- 
ment control, meant that sometimes the LVDTs measuring the displacement of the greater dis- 
placing top flange - see Figure [Ijb) - went out of range very quickly after the secondary instability 
was triggered. However, there were no such problems associated with the bottom flange measure- 
ments and so these were the primary values used for comparison purposes, since the theoretical 
model gives both Us and (j) directly. 



4 Numerical results and validation 

4.1 Cellular buckling 

The system of nonlinear ordinary differential equations p4|) - (|35p . subject to boundary conditions 
from equations p6p - ([55)) and integral equations from ([5^ - (|42p are solved using the well-known 
and tested numerical continuation package Auto [14]. For illustrative purposes. Figures [9HTT] 
present results from the variational model for test specimen 1. Figure [9] shows a plot of the 
(a) normalized moment ratio m, which is defined as the ratio M/Mcr, and (b) the normalized 
local buckling displacement amplitude (w{L/2)/tf) of the vulnerable part of the compression 
fiange versus the normalized lateral displacement of the bottom flange, {us — Ut^)lb. The graphs 
in (c) and (d) show the relationships between the "sway" and "tilt" components of the weak 
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Idealized beam in experiments: 



Schematic of experimental setup: 



^1 r- Lateral restraints — i ^\ 



"1 



l 



2^1 



XL 




Figure 7: Experimental rig designed to test fabricated I-beams under four-point bending. Although 
each beam was 4 m in length, the distance between the pinned supports was in fact 3.8 m. In each 
test, the lateral restraints were adjusted symmetrically to obtain the desired buckling length Le- 



axis centroidal displacement {qs and qt) and the torsional angle (q^ and g,-). A dotted line is 
superimposed on these graphs to show the Euler-Bernoulli assumption, where the sway and tilt 
amplitudes would be equal; this shows that the shear strains developed are small but not zero. 
Moreover, Figures IH^a) and (b) show a series of paths separated by a sequence of snap-backs 
with Figure [TU] presenting detailed graphs showing the corresponding numerical solutions beyond 
individual snap-backs for the local buckling functions. A distinctive pattern is clearly seen to 
emerge where the response passes from one path to the next, i.e. from CI to C2 to C3 and so on, 
in which each new path reveals a new local buckling displacement peak or trough. A selection of 
3-dimensional representations of the beams using the solutions for the paths CI, C3, C5 and C7 
are presented in Figure [TTl which include all components of LTB {us, (j), Utb and uu) and local 
buckling {w and u). 

As the response advances to path Cll, which has torsional rotations that are well beyond the 
scope of the model in terms of geometric considerations, the local buckle pattern is all but periodic 
and further loading would restabilize the system globally, assuming no permanent deformation has 
taken place. This global restabilization would occur as a result of the boundaries confining the 
spread of the buckling profile any further. Of course, if plasticity were present in the flange then 
any restabilization would be less significant and displacements would lock into plastic hinges. 

The phenomenon demonstrated currently and described above, where a sequence of snap-backs 
causes a progressive change from an initially localized post-buckling mode to periodic, has been 
termed in the literature as cellular buckling [15^ or snaking |16| . It has been found to be prevalent in 
systems where there is progressive destabilization and subsequent restabilization [26l [27] , such as in 
cylindrical shell buckling |T7| ITS' and kink banding in confined layers [T^, f^O] . In the fundamental 
studies concerning the model strut on a nonlinear foundation, the load oscillates about the Maxwell 
load, where the buckling modes progressively transform from localized (homoclinic) profiles to a 
periodic mode in a heteroclinic connection [351 US] • The oscillation in the strut model is attributed 
to the combination of nonlinearities in the foundation that have softening and hardening properties. 
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Array of LVDTs at beam midspan (cross-section view) 



Angle to account for 
beam vertical displacement 



LVDT: Top flange 



1 



LVDT: Vertical displacement 



Plate to account for 
beam lateral displacement 




^^ 



z. 



^^ 



Angle to account for 
beam vertical displacement 

LVDT: Bottom flange 

M 



Figure 8: Cross-section of the tested beam showing the relative positions of the LVDTs. Angles 
and plates were clamped as shown to both flanges such that the overall beam displacement still 
allowed the relative flange and vertical displacements to be evaluated. 

In the present context, as in the case for sandwich struts [151, the destabilization is derived from 
the interaction of instability modes with the restabilization arising from the inherent stretching 
that occurs during plate buckling due to large deflections, which accounts for its significant post- 
buckling stiffness |7]. Moreover, since the moment ratio m has a decreasing trend rather than 
oscillating about a fixed value, it is suggested that the destabilization is inherently more severe 
than the restabilization for the present case. 

4.2 Comparison with existing experiments 

Work conducted by Cherry ilOj , which focused on the overall buckling strength of beams under 
bending that had locally buckled flanges, presented a series of test results and proposed a theoret- 
ical estimate of the post-buckling strength. The theoretical approach was based on the effective 
width of the locally buckled flange which originated in [30]. However, the model presented in [10] 
was limited because of the assumption that both outstands of the compression flange behaved 
symmetrically. Nevertheless, the tests that were presented therein provide valuable data for the 
wavelengths of the local buckling mode in the compression flanges that were measured for four 
separate doubly-symmetric I-beams, with properties as presented in Table [31 the data are used 
for comparison purposes in the current study. Since the buckling mode predicted by the cur- 



Beam 


h (mm) 


b (mm) 


tuj (mm) 


tf (mm) 


E (kN/mm2) 


A 


72.3 


76.6 


3.68 


1.60 


66.47 


B 


74.0 


76.6 


3.78 


2.11 


64.53 


C 


72.2 


76.3 


5.13 


1.89 


64.43 


D 


71.5 


75.9 


4.65 


1.90 


66.19 



Table 3: Cross-section properties of beams tested and results presented in jlO^. Three different 
lengths of beams were tested under bending for each cross-section (A-D) ranging from 1 m to 2 m. 

rent analytical model is not necessarily periodic, but tends to approach this quality in the far 
post-buckling range, comparisons between the tests in [10 and the current model can be made 
when the profile for w exhibits periodicity throughout the beam length. Figure [T^] shows how 
the wavelength is defined from the post-buckling mode that has a central portion which is close 
to periodic. Table U shows the range of wavelengths obtained from the numerical solution of the 
system of equilibrium equations presented in ij2.2.6l The current model was run for a range of 
lengths between 1-2 m since this was the range for which the vast majority of tests presented 
in |10| were conducted. Apart from the beams with cross-section B, the comparisons are very 
encouraging; the discrepancies between the model and the tests are attributed to boundary effects 
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0.2 0.4 0.6 0.8 
(u —u )/h 




0.2 0.4 0.6 O.i 

{u —u )/b 




0.01 0.02 0.03 0.04 0.05 

9, 



2 








(d) 












Cll ■' 


1.5 










y^ 


1 








y^ 


- 


0.5 

fV 


y 


/- 


CAy/' 




^ 



0.5 1 

q (rad) 



1.5 



Figure 9: Numerical equilibrium paths for test specimen 1. Graphs of (a) normalized moment 
ratio m and (b) local buckling displacement amplitude Wmax/i/ versus lateral displacement of 
the bottom flange are shown. Points C and S show the critical and secondary bifurcation points 
respectively. Note that the sequence of paths, separated by snap-backs that mark the appearance 
of a new local buckling "cell", are denoted as "Cn" where n is the cell number. The graphs in (c) 
and (d) show the relationships between the generalized coordinates defining LTB during interactive 
buckling; note that the dotted lines represent the relationship between the respective quantities 
assuming Euler-Bernoulli bending, which show that the developed shear strains are small. 



affecting the results of the analytical model. It has been seen in the cellular buckling results earlier 
in this section (Figures IMTTj) that, as each buckling cell develops, the buckling "wavelength" A, 
see Figure [T2I drops until the buckling profile eventually tends to true periodicity and the moment 
M tends to a constant. For the numerical results from the current model that overestimated the 
wavelength, lack of convergence became an issue and the local buckling profile w was still showing 
remnants of the decaying tails near the boundaries, which are the signatures of homoclinic be- 
haviour. Hence, those particular comparisons are perhaps not entirely representative of the actual 
response predicted by the model. 



4.3 Results from current experiments 

For each of the physical experiments performed in the present study, see ^ testing proceeded 
to failure and all of the beams exhibited an unstable response once interactive buckling was 
triggered. A selection of photographs is presented in Figure [13] which show the beams from a 
variety of directions while they were undergoing interactive buckling. In tests 2 and 6, there 
was visual experimental evidence of cellular buckling; Figure [14] shows a sequence of photographs 
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Figure 10: Numerical solutions of the displacements of local buckling for the tip (x = —b/2) of the 
vulnerable flange of test specimen 1: out-of plane displacement w (left) and in-plane displacement 
u (right). Individual solutions on equilibrium paths CI to C7 (defined in Figure [9]) are shown in 
sequence from top to bottom respectively. All dimensions are in millimetres. 

before and after the principal instability showing a new local buckling peak appearing soon after 
the initial one. Table [S] presents the results and their comparison with the individual buckling 
modes. The maximum applied moment in the experiments A/max is presented as a ratio of the 
theoretical critical moment Af *-" calculated from the appropriate critical mode given in Table [2l 
whether LTB or local buckling. The local buckling profile was determined by marking (as seen 
in Figure [TSf d)) and measuring between adjacent peaks of the local buckling displacement over 
the length of the vulnerable part of the compression flange while the beam was still loaded but 
well after the peak moment had been applied. The interactive mode was clearly modulated in 
each case, with the peak amplitudes from local buckling decaying towards the lateral restraints; 
this was particularly notable in the cases where LTB was critical since the number of peaks was 
visibly fewer. In each test, two or three peaks of the local buckling mode exhibited significant 
plastic deformation almost immediately after the interactive mode was triggered; it was adjacent 



Beam 


Test wavelength [lOj 


Model wavelength A (mm) 


Error range 




(mm) 


Minimum 


Maximum 


(%) 


A 


126.2 


127.0 


133.4 


-hO.6-^+5.7 


B 


124.0 


153.0 


156.7 


-1-23.4^-^26.4 


C 


108.5 


103.1 


132.6 


-4.9^-1-22.3 


D 


109.7 


102.1 


141.2 


-7.0^-^28.7 



Table 4: Comparisons of buckling wavelengths between the tests in ^lOj and the current model. 
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(a) CI: m = 0.9998 



(b) C3: m = 0.9900 
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(c) C5: m = 0.9718 



X 

(d) C7: m = 0.9414 



Figure 11: Numerical solutions of the system of equilibrium equations visualized on a represen- 
tation of the actual beam used in test 1. All deformation components: w, u, Ug, 0, Utt and Utt 
are included with all dimensions in millimetres. The results are shown for individual points on 
paths CI, C3, C5 and C7 (Figure |9|) with the specific moment ratio m given. Note how the local 
buckling mode develops and how the "wavelength" of the local buckling profile within the central 
portion of the flange in more compression changes as more cells develop. 
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Figure 12: Definition of local buckling wavelength A from results for w from the variational model. 







Initial post-buckling 


Local buckling mode observations 


Test 


M,„ax/M^ 


drop in moment (%) 


No. of visible 
buckling peaks 


Extent 

(% of Le) 


1 


1.03 


15 


5 


25.4 


2 


0.90 


7 


7 


36.5 


3 


0.93 


12 


7 


40.6 


4 


1.07 


24 


8 


52.4 


5 


1.00 


14 


9 


65.8 


6 


1.05 


15 


9 


69.2 



Table 5: Results from the experimental programme in terms of the maximum moment and the 
local buckling profile. The final column is a measure of localized nature of the flange buckle - the 
smaller the number, the more it was localized. 



to these peaks where the buckling wavelength was, in general, measured to be the smallest values. 
Another notable feature shown in Table [5] is the immediate proportional drop in the moment 
once the interactive mode had been triggered. As would be expected from the literature [31], the 
largest drop occurred in test 4 where the critical modes had been practically simultaneous. It is 
also noteworthy that the tests with identical buckling lengths (tests 1 and 2) showed very different 
peak moments and moment drops. A rational hypothesis can be devised for this by postulating 
that the beam in test 2 contained more geometric imperfections than the beam in test 1. This 
would not only account for the smaller maximum moment measured in test 2, but also for its 
smaller relative moment drop and its lower residual moment in the post-buckling range [9]. 

4.3.1 Comparisons with variational model and discussion 

Figure [15] presents normalized plots of the applied moment m versus the measured and normalized 
lateral displacement of the bottom flange, {ug — u^)/b. Test 4 gives clearly the best comparison 
in terms of the correlation between the post-buckling response of the actual beam and the model 
prediction. Tests 1 and 2 also show good basic agreement with the theory; test 1 showing that the 
post-buckling unloading resembles the theory quite well, while test 2 shows that the instability 
is triggered at a similar value of the lateral displacement predicted by the theory (see Table [6]) . 
Tests 5 and 6 clearly peak at or marginally above the local buckling critical moment, as predicted 



Test 


[Us - 


- u„)/b 


Local buckling wavelengths A (mm) 




Expt 


Theory 


Expt range 


Expt average 


Theory (minimum) 


1 


0.169 


0.134 


185 -^ 221 


203 


200 


2 


0.146 


0.134 


171 ^ 250 


195 


200 


3 


0.061 


0.110 


161 -^ 242 


203 


150 


4 


0.068 


0.083 


161 -^ 249 


206 


153 


5 


0.066 


0.052 


167 ^ 247 


206 


163 


6 


0.061 


0.058 


158 -^ 225 


195 


240 



Table 6: Comparisons of the experimental results with the variational model in terms of the values 
of the bottom flange displacement at the peak moment point and the local buckling wavelengths. 

from linear analysis. However, in a similar way to test 2, the instability is triggered at a lateral 
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(a) Test 1: Close up of local buckling waves 



(b) Test 3 




(c) Test 6 



(d) Post-testing beams: Tests 1, 3, 4 and 5 



Figure 13: Selection of photographs from the experimental programme: (a)-(c) all showing inter- 
active buckling and (d) shows four of the beams and their locally buckled flanges that show plastic 
deformation. 



displacement that correlates well with the prediction from the variational model. For test 6, the 
variational model yields a lower critical moment than the Mcr value for LTB, which triggers a 
quasi-local buckling mode. However, as stated earlier, a distinct and accurate local buckling mode 
can only be modelled with additional displacement functions in the current framework so this 
particular result needs to be interpreted with some caution. Test 3 could be considered to be 
an outlier, but the measured response would imply, in a similar way that was discussed above 
regarding test 2, that the level of geometric imperfections in this beam was higher than the other 
tested beams (1, 4, 5, 6). Hence, the measured instability moment is less, the unloading proportion 
is less and the response is practically parallel to the model curve, which in fact is encouraging. 

In terms of the local buckling wavelengths, these are compared to the wavelength of the buckling 
profile obtained from the variational model as described in §4.21 Even though the theoretical 
results seemed to be influenced by effects close to the boundary (particularly in test 6), hence the 
variability in the predictions, the general correlation between the experiments and theory is good. 
The apparent confirmation that the post-buckling behaviour of an I-beam under pure bending is 
cellular when global and local instability modes interact nonlinearly poses the following question: 
is this phenomenon prevalent in other thin-walled structural components that are known to suffer 
from overall and local mode interaction? Compressed stringer-stiffened plates |7] and I-section 
struts [51 are prime examples of other components where local and global mode interactions are 
known to occur. Further research is obviously required to determine the answer. 



Article submitted to The Royal Society 



TeX Paper 



Cellular buckling in I-beams 



19 










(a) Test 2 (initial) 
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(b) Test 2 (LTB and 1st local buckle formation) 












(c) Test 2 (2nd local buckle formation) 
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(d) Test 6 (initial) 
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(e) Test 6 (LTB and 1 st local buckle formation) 
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(f) Test 6 (2nd local buckle formation) 

Figure 14: Evidence of cellular buckling. Two sequences of three photographs are shown of tests 
2 (a)-(c) and 6 (d)-(f) respectively. Photographs (a) and (d) respectively show the pre-buckling 
state; (b) and (e) respectively show the initial post-buckling with one significant peak at midspan; 
(d) and (f) respectively show a newly developed local buckling peak in the top flange. 
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Figure 15: Moment ratio m versus bottom flange lateral displacement {us — Uw)/h for all the exper- 
iments with the variational model predictions superimposed, denoted as "Theory" and evaluated 
by Auto. Points C and S refer to the critical and secondary instability points obtained from the 
model. Note that local buckling was the theoretical critical instability mode for Tests 4-6. 
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5 Concluding remarks 

The current work identifies an interactive form of buckling for an I-beam under uniform bending 
which couples a global instability with local buckling in one-half of the compression flange. In 
contrast to earlier, more numerical, work |13 1 124 ] . cellular buckling, the transformation from a 
localized to an effectively periodic mode, is predicted theoretically for the purely elastic case and 
evidenced in physical tests. The model compares well both qualitatively and quantitatively with 
the observed collapse of a beam that undergoes the interaction under discussion that involves 
global, local, localized and cellular buckling. The localized buckle pattern first appears at a 
secondary bifurcation point which immediately destabilizes a portion of the compression flange; 
as the deformation grows, the buckle tends to spread in cells until eventually it restabilizes when 
the localized buckling pattern has become periodic after a sequence of snap-backs. 

Experimentally, the process is unstable and so this sequence occurs rapidly even under rigid 
loading with the local buckling cells being triggered dynamically. This highlights the practical 
dangers of the modelled and observed phenomenon; the interaction reduces the load carrying 
capacity, it therefore introduces an imperfection sensitivity that would need to be quantified such 
that robust design rules can be developed to mitigate against such hazardous structural behaviour. 
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